Phase diagram of Janus Particles 
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We deeply investigate a simple model representative of the recently synthesized Janus particles, 
i.e. colloidal spherical particles whose surface is divided into two areas of different chemical composi- 
tion. When the two surfaces are solvophilic and solvophobic, these particles constitute the simplest 
example of surfactants. The phase diagram includes a colloidal-poor (gas) colloidal-rich (liquid) 
de-mixing region, which is progressively suppressed by the insurgence of micelles, providing the first 
model where micellization and phase-separation are simultaneously observed. The coexistence curve 
is found to be negatively sloped in the temperature-pressure plane, suggesting that Janus particles 
can provide a colloidal system with anomalous thermodynamic behavior. 



An unprecedented development in particle synthesis is 
generating a whole new set of colloidal particles, char- 
acterized by different patterns, particle patchiness and 
functionalities [1. This extensive effort is driven by 
the attempt to gain control over the three-dimensional 
organization of self-assembled materials, characterized 
by crystalline, or even disordered, colloidal structures 
with desired properties. The recently synthesized Janus 
colloids [2j IH 0J [5j [H [7] — characterized by a surface di- 
vided into two areas of different chemical composition 
— offer a fine example of such versatility and can be 
reckoned as the paradigm of a system formed by incom- 
patible elements within the same unit structure which 
spontaneously drive a self-assembly aggregation process 
into complex superstructures. Not only this is ubiquitous 
in biological and chemical systems, but, when carried out 
into a controlled manner, it constitutes one of the most 
promising bottom-up process for the design of future ma- 
terials. Indeed the assembly behavior of Janus particles 
is receiving a considerable attention [HI l9l [TO] . 

With a proper choice of the chemico-physical surface 
properties, Janus particles can provide the most ele- 
mentary and geometrically simple example of surfactant 
particles [11], in which solvophilic and solvophobic areas 
reside on different parts of the surface of the same parti- 
cle, in a controlled analogy to what is often found in the 
protein realm, e.g. hvdrophobin[T2] and casein [T3]. 

The phase behavior of typical surfactant molecules 
(as opposed to colloids) has been studied in details in 
the past, due to their role in industrial processes and 
products [14]. It is well appreciated that, with the right 
balance of solvophilic and solvophobic moieties, surfac- 
tants can spontaneously self assembly into a large vari- 
ety of structures, including micelles, vesicles, lamellae. 
Often, an attractive interaction between the aggregates 
originates a macroscopic phase separation into a dense 
and a dilute micellar phase [15]. In rare cases, the sur- 
factant rich phase is characterized by micelles and the 
surfactant poor phase by unimers. One interesting ques- 
tion concerns the relation (or competition) between the 
gas-liquid phase separation and the process of micelles 
formation. While it is recognized that micelles formation 



can be seen as a sort of phase separation (in such a case 
the liquid is interpreted as an infinite size aggregate) it is 
rather unclear how the two phenomena really interrelate 
and if the same system can show a clear self assembly 
process in a certain region of temperatures T and densi- 
ties p and a clear gas-liquid phase separation somewhere 
else. Previous extensive search for models able to display 
this mixed behavior — carried out by Panagiotopoulos 
and coworkers [T6 ] [T7 | [T8] — has been unsuccessful, thus 
hampering a deep understanding of the interplay between 
aggregation and phase separation. In all studied cases, 
on tuning some model parameters, macrophase separa- 
tion becomes replaced by self-organization into isolated 
or interconnected micelles [T9] . 

In this Letter, we show that simply but realistically 
modeled Janus particles are ideal candidates for study- 
ing, within the same system, the interplay between aggre- 
gation and phase separation. Indeed, the system shows a 
clear phase-separation, characterized by the presence of 
a gas-liquid critical point and a micellar phase at lower 
temperatures. Micelles develop in the gas colloidal phase 
progressively acquiring stability, suppressing more and 
more the two-phase coexistence. Differently from all pre- 
viously studied models, this micelles rich gas phase coex- 
ists with a dense fully connected liquid phase. 

The Janus particles are described via the one-patch 
model introduced by Kern and Frenkel [20] with cover- 
age 0.5. In this effective model the two hemispheres of 
a hard-core particle of diameter a are considered respec- 
tively repulsive and attractive. When the segment con- 
necting the centers of two particles crosses the attractive 
hemisphere of both particles, the two particles then in- 
teract via a square- well potential with range A = 0.5a 
and depth uq. In all other orientations the interaction 
is hard-core. This model is inspired by the experimental 
system investigated in Ref . [6] , where the repulsive inter- 
action has an electrostatic origin and the attractive part 
is hydrophobic. In the following, a provides the unit of 
length and Uq the unit of energy. T is also expressed in 
unit of uq (i.e. Boltzman constant ks = 1). 

Results reported in this Letter are based on an ex- 
tensive numerical study of the model in a wide range 
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of T and p. We employ the most accurate numeri- 
cal methodologies to provide a detailed quantification of 
the phase behavior, including size-effect studies. Grand- 
canonical Monte Carlo (GCMC) simulations (for box size 
L = 10, 12, 15, 24a) have been implemented to locate the 
critical point and the relation between density and activ- 
ity, Gibbs ensemble simulations [21] to evaluate the co- 
existing phases and standard Monte Carlo (MC) simu- 
lations in the NVT and NPT ensembles for large sys- 
tem sizes (number of particle N = 5000) to elucidate the 
structure of the phases and the equation of state. Since 
self-assembly requires large attraction strengths (com- 
pared to ksT) and large systems (indeed aggregates can 
grow to very large sizes) accurate numerical studies are 
rather difficult. We have explicitly checked that the stud- 
ied system sizes are sufficiently large to avoid size effects. 
Large system sizes (N > 1000) have been studied also in 
GCMC and GEMC, especially at small densities to make 
sure that the simulation box contains several aggregates 
simultaneously. In all cases, translational and rotational 
moves consisting of a maximum random translation of 
±0.1<j and a maximum random rotation of ±0.1 rad of a 
randomly selected particle, have been implemented. De- 
pending on the MC method, insertion and deletion moves 
(or swap moves) have been attempted, on average, ev- 
ery 500 displacement moves while volume change moves 
(with maximum volume changes of the order of 0.5cr 3 ) ev- 
ery 100 translational moves. Extremely long simulation 
runs (of the order of 10 8 MC steps) have been performed 
to reach equilibrium. GEMC and GCMC simulations 
have been particularly painstaking, due to the difficulty 
of inserting/deleting, by a sequence of individual moves, 
large aggregates, requiring — for each studied low T state 
point — several months of computer time. 
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FIG. 1: Average potential energy (E) per particle along sev- 
eral isotherms vs. number density p = N/L 3 . Note that 
—2(E)/Nuo coincides with the average number of bonded in- 
teractions per particle. The thick lines connect the two coex- 
isting densities (full symbols) calculated via GEMC. 

Fig. [I] shows the average potential energy (E) per par- 



ticle along several isotherms. Below ksT/uo = 0.3, the 
density dependence of (E) becomes non monotonic (rem- 
iniscent of the shape of a van der Waals loop), already 
suggesting the possibility of a phase separation. Sev- 
eral features significantly differ from the standard gas- 
liquid behavior: the most striking difference is the fact 
that, close to the loop, particles in the low density states 
are characterized by E values significantly smaller than 
the one in the dense phase. At odd with the usual gas- 
liquid behavior, the position of the loop shifts to larger 
p and its width appears to shrink on cooling. A glance 
to the structure of the system across density ( Fig. [5| 
along different isotherms immediately clarifies the rea- 
sons for such behavior. While at higher T one observes 
the behavior typical of a fluid close to the gas-liquid crit- 
ical point, i.e. a spatially inhomogeneous density (or 
equivalently particles associated in highly poly-disperse 
clusters), at lower T the system self-assembles into well- 
defined, essentially spherical, clusters which appear to 
maintain their character in a wide region of p. The most 
frequent clusters are composed by single-layer micelles 
(in which about ten particles orient their attractive part 
inside the aggregate, Fig. [3|(a) ) and double-layer mi- 
celles (or vesicles, in which about 40 particles organize 
into an inner and an outer shell forming a bilayer which 
exposes to the outside only the hard-core part, prevent- 
ing any attraction between distinct aggregates, Fig. [3} 
(b)). The double-layer structure becomes the dominant 
one on cooling and/or on increasing p. The organization 
of the system into these structures has the net effect of 
breaking the thermal correlations which propagate via 
inter-particle attraction, effectively preventing the devel- 
opment of long range critical fluctuation. This effect is 
evident at p = 0.1 (Fig. J2^c),(h),(m)), where a percolat- 
ing network of bonded particles is found only at T = 0.30. 
The insurgence of an energy driven formation of micelles 
effectively destabilizes the demixing of the system into a 
colloidal-poor (gas) phase coexisting with a colloidal-rich 
(liquid) phase. 

To provide evidence of the existence of a critical point, 
of a truly phase-coexistence and to evaluate the appropri- 
ate phase boundaries , we have performed GCMC sim- 
ulations to locate the T and chemical potential values 
for which the density fluctuations behave as the mag- 
netization in the Ising model [22 (providing the location 
of the critical point) and GEMC simulations, to calcu- 
late the densities of the coexisting phases. Results are 
reported in Fig. |1J Differently from the standard behav- 
ior, the region of phase coexistence which usually widens 
on cooling, here shrinks and shifts to the right, leaving 
a wide p window where a homogeneous gas-like phase 
composed by large aggregates [Fig. [2] (g) -(h); (m)-(o)] is 
thermodynamically stable. Micelles coexist with a liquid- 
like non-micellar phase [Fig. [5] (j) and (p)], supporting 
the evidence that the phase-separation is not related to 
a critical phenomenon associated with inter-micelles in- 
teractions. The densities of both coexisting phases in- 
crease on cooling, progressively approaching each other, 
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FIG. 2: Snapshot of an equilibrium typical configuration for a 5000 particles system at ksT/uo — 0.3 (a-e), ksT/uo = 0.27 
(f-j) and ksT/uo = 0.25 (k-q) for selected densities. The green-coded area is attractive (square- well) , the red-coded is repulsive 
(hard-sphere). The scale changes according to the density. 





FIG. 3: Typical shape of the aggregates of size 10 (a, micelles) 
and 40 (b, vesicles). 



suggesting the possibility of a lower T consolute point. 
Unfortunately, both the gas of micelles and the liquid 
phase become metastable in respect to the formation of 
a ordered (crystalline, lamellar) phase; the onset of crys- 
tallization, on the time scale probed by the numerical cal- 
culations, pre-empts the possibility of exploring such pos- 
sibility. A snapshot of the system ordered in the lamellar 
phase is reported in Fig. |2|q). 

The energetic and entropic balance associated with the 



transition from the micelles gas to the liquid phase reveals 
some important information. As previously alluded while 
discussing Fig. [T] differently from the usual gas-liquid 
behavior, (E) is higher in the liquid phase than in the 
micellar gas phase. Hence, despite the fact that the gas 
phase is stabilized by the translational entropy of the mi- 
celles, at coexistence the liquid phase is probably more 
(orient ational) disordered than the gas one. Insights on 
the larger orientational order in the gas phase as com- 
pared to the liquid one, can be obtained by computing 
the distribution of relative orientations between all pairs 
of bonded particles. More precisely, for all bonded pairs, 
we have evaluated the distribution of the scalar product 
ni • n 2 where ni and n 2 are the two unit vectors indi- 
cating the location of the patch center in each particle 
frame (see arrows in Fig. 5|. The distribution P(ni • n 2 ) 
is expected to show well defined peaks for an ordered 
state and to be flat in a completely disordered state. As 
shown in Fig. [5j P(ni • n 2 ) is indeed significantly more 
structured in the micelles-rich gas phase with respect to 
the liquid phase. 

To provide a definitive confirmation of the lower en- 
tropy of the gas phase (S g ) as compared to the liquid 
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FIG. 4: Phase diagram of the Janus model in the T — p (a) 
and P — T (b) planes. Symbols are simulation results, lines 
are guide to the eyes. 
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FIG. 5: Distribution of the scalar product ru • n2 evaluated 
over all pairs of bonded particles in a gas (pa 3 = 0.65) and 
in a liquid (pa 3 = 0.675) state point, both at kT/uo = 0.25. 
Cartoons provide a sketch of the relative orientations of the 
bonded particles. Arrows indicate the corresponding unit vec- 
tors. 



one (Si), we evaluated the equation of state at several 
T by mean of MC simulations in the PNT ensemble. 



These results, coupled with the coexisting densities eval- 
uated via GEMC, provide a method to calculate the 
slope (dP/dT\ x ) of the coexistence curve in the P — T 
plane. By virtue of the Clapeyron equation, dP/dT\ x is 
a measure of AS/ AV , where AS and AV are the differ- 
ences in entropy and volume of the two coexisting phases. 
Since V g > V\ a positive (negative) slope implies that 
S g > (<)Si. Fig. [41(b) shows that the coexistence curve 
in the (T — P) plane has an anomalous behavior. Dif- 
ferently from the usual gas-liquid coexistence, here the 
curve is negatively sloped, confirming that S g < Si, de- 
spite the significantly larger translational contribution. 
The explanation for this odd behavior can be traced back 
to the orientational ordering of the particles in the gas 
micellar phase and to the disordered orientations which 
characterize the liquid phase (Fig.|5|. As a result, the ori- 
entational disorder gain in the fluid phase compensates 
the energetic penalty. The negatively sloped gas-liquid 
coexistence curve indicates also that, on cooling along 
a constant (osmotic) pressure path, the density of the 
system decreases at the transition. 

The results presented here have far-reaching conse- 
quences in many respects. From the theoretical side, our 
work constitutes the first study displaying a clear gas- 
liquid critical point which evolves into phase coexistence 
between a cluster and a liquid phase through a compe- 
tition between phase-separation and micelles formation. 
In addition, the system studied here behaves as anoma- 
lous substance in which the dense phase is more disor- 
dered than the low-density phase and which, in analogy 
with the negatively sloped liquid water- hexagonal ice co- 
existence curve, expands on cooling along isobars. The 
simplicity of the model and its rich and unconventional 
phase diagram begs for experimental verifications. De- 
spite the investigated model has a range (0.5a) larger 
than the interaction ranges typical of /im-sized colloids 
(usually of the order of 0.05-0.2 a), ranges comparable 
to 0.5<j can be realized with nanoparticles or building on 
the recently explored Casimir critical forces [23], [24] . 
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